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Abstract 

We compute the quark-mass dependence of the baryon octet and decuplet masses 
within the X"MS scheme at the one-loop level. The results are confronted with 
recent lattice QCD simulations of the MILC collaboration. A fair reproduction of 
the quark-mass dependence as suggested by the lattice simulations is obtained up 
to pion masses of about 700 MeV. Based on the chiral one-loop results we predict 
that the dependence of the baryon octet and decuplet masses on the quark-masses is 
discontinuous. Typically the pion-mass dependence is smooth up to about 300 MeV 
only. This is a consequence of self consistency imposed on the partial summation, i.e. 
the masses used in the loop functions are identical to those obtained from the baryon 
self energies. The 'mysterious' quark-mass dependence of the H mass predicted by 
the MILC collaboration is recovered in terms of a discontinuous chiral extrapolation. 



1 Introduction 



The present-day interpretation of QCD lattice simulations requires a profound 
understanding of the dependence of observable quantities on the light quark 
masses. A powerful tool to derive such dependencies is the chiral Lagrangian, 
an effective field theory based on the chiral properties of QCD. The applica- 
tion of strict chiral perturbation theory to the SU(3) flavor sector of QCD 
is plagued by poor convergence properties for processes involving baryons 
[1,2,3,4,5]. Thus it is important to establish partial summation schemes that 
enjoy improved convergence properties and that are better suited for chiral ex- 
trapolations of lattice simulations. It was suggested by Donoghue and Holstein 
[6,7,8] to improve the convergence properties by introducing a finite cutoff into 
the heavy-baryon effective field theory of Jenkins and Manohar [9]. Another 
scheme was suggested in [10,11,12], the construction of which was guided by 
keeping manifestly covariance, analyticity and causality. The computation of 



Preprint submitted to Elsevier Science 



9 February 2008 



the baryon octet and decuplet masses in this scheme was worked out recently 
indicating convincing convergence properties of the chiral loop expansion [12]. 

For studies of chiral extrapolations of the nucleon mass within the chiral SU(2) 
Lagrangian we refer to the work by Procura, Hemmert and Weise [16]. The 
latter study applied the partial summation scheme of Becher and Leutwyler 
[17]. It was based on simulations of the CP-PACS collaboration that used dy- 
namical u- and d-quarks only [18]. The application of the scheme by Becher 
and Leutwyler [17] to the chiral SU(3) Lagrangian is questioned by poor con- 
vergence properties as was demonstrated by EUis and Torikoshi at hand of the 
baryon octet masses [4]. The recent work of Pascalutsa and Vanderhagen [19] 
suggested to take serious the partial summation implied by the extended on- 
mass shell scheme (EOMS) introduced by Gcgclia and Japaridze [20]. Their 
computation addressed so far the chiral extrapolation of the nucleon and iso- 
bar mass only. 

It is the purpose of the present work to confront the results of [12] to recent 
lattice QCD simulation of the MILC collaboration [13,14], that use dynam- 
ical u-,d- and s-quarks in the staggered approximation. We do not aim at 
a quantitative extrapolation of the lattice simulation, which would require a 
continuum limit extrapolation and a quantitative control of systematic errors. 
Rather, we want to use the available partial lattice results to make predictions 
for future QCD lattice simulations and possibly obtain rough constraints on 
some chiral parameters. An analysis of the recent MILC results, similar in 
spirit, was undertaken by Prink and MeiBner [15] based on the cutoff scheme 
of Donoghue and Holstein. 

Adjusting the values of the counter terms to the physical masses we predict 
a discontinuous dependence of the baryon masses on the pion mass. This is a 
consequence of self consistency imposed on the partial summation approach, 
i.e. the masses used in the loop function are identical to those obtained from 
the baryon self energy. The latter is a crucial requirement since the loop func- 
tions depend sensitively on the precise values of the baryon masses. Our results 
may explain the mysterious quark-mass dependence observed by the MILC 
collaboration for the S mass [13,14]. 



2 Chiral interaction terms 



We collect the terms of the chiral Lagrangian that determine the leading orders 
of baryon octet and decuplet self energies [21,22]. Up to chiral order Q'^ the 
baryon propagators follow from 
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We assume perfect isospin symmetry through out this work. The fields are 
decomposed into isospin multiplets 



(2) 



with the Gell-Mann matrices, Aj, and the isospin doublet fields K — {K^, K^Y 
and S = The isospin Pauli matrices a = (cti, (T2, Cs) act exclusively 

in the space of isospin doublet fields {K, N, 5) and the matrix valued isospin 
doublet a, 



a' ^ -^{X^ + i A5, Ae + i A7) , r = (Ai, A2, A3) . 



(3) 



The evaluation of the baryon self energies to order probes the meson-baryon 
vertices 
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where we apply the notations of [11]. We use / = 92.4 MeV in this work. 

The values of the coupling constants F, D, C and H may be correlated by a 
large- Ai'c operator analysis [23,24,25]. At leading order the coupling constants 
can be expressed in terms of F and D only. We employ the values for F and 
D as suggested in [26,11]. All together we use 



0.45, D = 0.80, H^9F-3D, C^2D, 



(5) 
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in this work. We take the parameter Z = 0.72 from a detailed coupled-channel 
study of meson-baryon scattering that was based on the chiral SU(3) La- 
grangian [11]. The latter parameter is an observable quantity within the chi- 
ral SU(3) approach: it contributes at order to the meson-baryon scattering 
amplitudes but cannot be absorbed into the available counter terms [11]. 



3 Chiral loop expansion at the one-loop order 

We briefly recall the results of [12]. A summation approach was defined by 
performing a chiral loop expansion rather than a strict chiral expansion: for a 
given truncation of the relativistic chiral Lagrangian we take the loop expan- 
sion that is defined in terms of the approximated Lagrangian seriously. The 
number of loops we would consider is correlated with the chiral order to which 
the Lagrangian is constructed. A renormalization based on the Passarino- 
Veltman reduction was devised that installs the correct minimal chiral power 
of a given loop function. The residual dependence on the renormalization scales 
is used to monitor the convergence properties of the expansion and therewith 
estimate the error encountered at a given truncation. 

The loop contribution to the baryon octet and decuplet masses read: 
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The sums in (6, 7) extend over the intermediate Goldstone bosons {Q e [8]) 
baryon octet {R G [8]) and dccuplet states {R G [10]). The various couphng 

( B) 

constants Gq^; arc determined by the parameters F, D.C, H. They are hsted 
in [12]. We emphasize that (6, 7) depend on the physical meson and baryon 
masses mg and Mr. This defines a self consistent summation since the masses 
of the intermediate baryon states in (6, 7) should match the total masses. At 
order the latter are the sum of the tree-level contributions linear in the 
parameter bo,b]:),bF or do, do and the loop contribution (6, 7). 

The mesonic tadpole Iq enjoys a logarithmic dependence on the ultravio- 
let renormalization scale fiuv and the one-loop master function Iqr a linear 
dependence on the infrared renormalization scale fijR. By construction the 
results (6, 7) arc necessarily consistent with all chiral Ward identities as dis- 
cussed in [12]. As emphasized by the authors the parameters ho, hi), hp and 
do, do are strongly dependent on the infrared scale ^jr. Applying a further 
chiral expansion to (6, 7) it was demonstrated that the physical masses are 
scale independent as they should be. However, the convergence properties re- 
flect the choice of the renormalization scales. Good convergence properties can 
only be expected for natural values thereof. For the infrared scale hir ~ Q a 
natural window 350 MeV < jiiR < 550 MeV was suggested in [12]. Keeping 
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the partial summation as defined by (6, 7) a residual dependence on the renor- 
malization scales remains. This is analogous to the residual cutoff dependence 
of the scheme of Donoghue and Holstein [6] . As long as such dependencies are 
small and decreasing as higher order terms are included they do not pose a 
problem, rather, they offer a convenient way to estimate the error encountered 
at a given truncation. 



4 Qu£trk-mass dependence of the baryon masses 

We discuss the implications of (6, 7) for chiral extrapolations of lattice sim- 
ulations of the baryon masses [13,14]. Since we are not aiming at a chiral 
extrapolation of the lattice simulations down to physical quark masses we ad- 
just part of the parameters to empirical data directly. The goal of the present 
study is a qualitative understanding of the quark-mass dependence of the 
baryon masses. 

For given values of the infrared and ultraviolet renormalization scales the pa- 
rameters ho, hp and do are fitted to the mass differences of the octet states and 
decuplet states. The absolute mass scale of the octet and decuplet states can 
be reproduced by appropriate values of the bare baryon masses. This proce- 
dure leaves undetermined the two parameters and d^. Good representations 
of the physical baryon masses can be obtained for a wide range of the latter 
parameters. The parameter 6o may be used to dial a given pion-nucleon sigma 
term at physical pion masses. Similarly the unknown parameter do m^-y be 
determined as to reproduce a given pion-delta sigma term. 

In this work 6o and d^ are adjusted as to reproduce the baryon octet and 
decuplet masses in the SU(3) hmit at ~ 690 MeV. The MILC simulations 
suggest the values Mfgj ~ 1575 MeV and M[io] ~ 1710 MeV. This procedure 
is biased to the extent that it assumes that the chiral one-loop results will be 
applicable at such high quarks masses. On the other hand as long as there 
are no continuum limit results of the MILC collaboration available this is an 
economical way to minimize the influence of lattice size effects. The latter are 
expected to be smaller at large quark masses. The procedure may be justified in 
retrospect if it turns out that the extrapolation recovers the behavior predicted 
by the lattice simulation. 

The parameters used in this work are collected in Tab. 1 together with the 
implied masses of the baryon octet and decuplet states. A fair representation 

of the physical baryon masses is obtained. For simplicity the values quoted in 
Tab. 1 correspond to a computation where the intermediate baryon masses 
are put to their empirical values. Since the resulting masses are very close 
to the physical masses this is well justified. As discussed in detail in [12] the 
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a^N [MeV] 


52.4 
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a^-p [MeV] 


384.0 
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UK-n [MeV] 


359.4 


354.8 


361.5 


Table 1 



The parameters are fitted so that (6, 7) reproduces the baryon masses at physical 
pion masses as well as the SU(3) limit values Mjg] ~ 1575 MeV and Mjio] — 1710 
MeV at ~ 690 MeV. We use i^uv = 800 MeV. The masses are decomposed into 
their chiral moments. 



parameters bQ^^ p and (io,_D show a strong dependence on the infrared scale 
HiR. In contrast the physical baryon masses suffer from a weak dependence 
only. For natural values of the infrared scale the chiral expansion appears well 
converging as indicated by the decomposition of the baryon masses into their 
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Fig. 1. The pion mass dependence of baryon octet and decuplet masses predicted 
by the chiral loop expansion taking the parameters of Tab. 1. The hnes represent 
the masses for the infrared scale put at fijR = 450 MeV. The solid squares are the 
simulation points of the MILC collaboration. 



moments. 



We emphasize that the values of the pion-nucleon sigma term shown in Tab. 1 
are an immediate consequence of the MILC simulation of the SU(3) symmetric 
point defined by ~ 690 MeV as discussed above. No further results of the 
MILC simulation are used. Taking the residual scale dependence of the pion- 
nucleon sigma term as a naive error estimate we obtain a-jrN = 54±2 MeV. The 
latter value is somewhat larger than the canonical value of Gasser, Leutwyler 
and Sainio cr-^N — 45 MeV [27]. Our value is, however, quite compatible with 
a recent analyses of Frink and Meissner who suggest aj^N — 52 MeV [15] . Also 
Pascalutsa and Vanderhagen [19] arrive at a somewhat larger value a-j^N — 57 
MeV based on a chiral SU(2) analysis of the MILC simulation points. 

We turn to the pion-mass dependence of the baryon octet and decuplet masses. 
Our results are collected in Fig. 1 based on the parameters of Tab. 1. We show 
results for the particular choice fijR = 450 MeV. It should be emphasized that 
the baryon masses are a solution of a set of coupled and non-linear equations 
in the present scheme. This is a consequence of self consistency imposed on the 



partial summation approach. The latter is a crucial requirement since the loop 
functions depend sensitively on the precise values of the baryon masses. As a 
consequence of the non-linearity for a given parameter set there is no guaran- 
tee for a unique solution to exist, nor that solutions found are continuous in 
the quark masses. Indeed as illustrated by Fig. 1 the pion-mass dependence 
predicted by the chiral loop expansion is quite non-trivial exhibiting various 
discontinuities. Upon inspecting the energy dependence of the baryon self en- 
ergies this is readily understood: the baryon self energies, as implied by (6, 7), 
allow in certain cases for multiple solutions of the mass equation. 



^^M+Eb{^), (9) 

even though the self consistency requirement singles out a particular solution 
at 1^ = Mb- Such a scenario appears inconsistent at first: the meson baryon 
loop functions have to be evaluated with all states provided by (9). In contrast 
the computation of the meson-baryon loop functions (6, 7) takes the existence 
of uniquely defined baryon states for granted. In addition it is assumed that 
the mass poles of the baryon propagators have residuum one. Thus our com- 
putation considers implicitly additional counter terms of the form: 



^wave- function = Cb (j^ - M^) , ( 10) 

where Mb denotes the baryon masses that solve the self consistent system. If 
the term (10) is added to the self energies it does not affect the solutions of 
the self consistent system. However, it does affect the energy dependence of 
the final baryon self energy. By adjusting the parameters Cb, it is possible to 
arrive at uniquely defined baryon states. 

Within the scheme described above it is intuitive that for a given set of pa- 
rameters multiple solutions of the self consistent system may exist. A priori it 
is not evident which of the possible solutions one should select. Our strategy 
is to continuously evolve the solution established at physical pion masses to 
larger pion masses. Similarly we evolve down to smaller pion masses the so- 
lution that reproduces the lattice simulations at the SU(3) symmetric point 
defined by rriT^ ~ 690 MeV. Both solutions reach end points beyond which we 
do not find any continuation. In the region of intermediate pion masses we, 
however, find another solution, suggesting a discontinuous dependence of the 
baryon masses on the quark masses of QCD. One may expect that the physical 
solutions are those of lowest mass always. The three branches discussed above 
are shown in Fig. 1 and confronted with the simulation points of the MILC 
collaboration [13,14]. 

It is striking to see that we reproduce the 'mysterious' pion-mass dependence 
of the S mass, i.e. the quite fiat behavior which does not seem to smoothly 
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approach the physical mass. Given the present uncertainties from finite lat- 
tice spacing, the staggered approximation and the theoretical uncertainties 
implied by higher order contributions, we would argue that we arrive at a fair 
representation of the lattice simulation points for all baryons with some reser- 
vation concerning the nucleon. We note that the parameter set used before in 
[12] leads to results qualitative similar to those shown in Fig. 1 for the baryon 
octet masses. While the description of the nucleon is improved showing no 
discontinuity at large pion masses, the sizeable jump of the S mass at small 
pion masses remains. However, the MILC simulation points of the isobar mass 
disfavor the latter choice. An overestimate of the lattice simulations of the 
isobar mass at larger quark masses would be the consequence. 

Incorporating the many counter terms offered by the chiral Lagrangian it is 
reasonable to expect that the latter will further improve the picture. However, 
as long as there is no detailed analysis available that performs the continuum 
limit and estimates the uncertainty from the staggered fermion approximation 
there is not much point considering the counter terms. 



5 Summary 

We evaluated the baryon octet and decuplet self energies at the one-loop level 
applying the relativistic chiral Lagrangian. Adjusting the counter terms an 
excellent representation of the physical baryon octet and decuplet masses is 
obtained. The size of the pion-nucleon sigma term was estimated by including 
in the fit the baryon octet and decuplet masses at the SU(3) symmetric point 
defined by ~ 690 MeV. For the latter the MILC collaboration suggests 
the values M^gj ~ 1575 MeV and M[io] ~ 1710 MeV. As a result we predict 
Ut^n — 54 MeV. The kaon-nucleon sigma terms are Ox-p — 380 MeV and 
(TK-n - 355 MeV. 

Given this scenario the pion-mass dependence of the baryon octet and de- 
cuplet masses were evaluated. The latter are a solution of a set of coupled 
and non-linear algebraic equations. This is a direct consequence of self con- 
sistency imposed on the partial summation, i.e. the masses used in the loop 
functions are identical to those obtained from the baryon self energies. As 
a striking consequence we predict a discontinuous dependence of the baryon 
masses on the pion mass. Typically the baryon masses jump at pion masses 
as low as 300 MeV. Most spectacular is the behavior of the S mass. At small 
pion masses it decreases with increasing pion masses. At a critical pion mass 
of about 300-400 MeV it jumps up to a value amazingly close to the predic- 
tion of the MILG collaboration. For all baryon masses for which the MILG 
collaboration published simulation points our results are reasonably close to 
the lattice estimates, given the present uncertainties from finite lattice spac- 
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ing, the staggered approximation and the theoretical uncertainties imphed by 
higher order contributions to the baryon self energies. 
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